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I report on a recent determination by the ALPHA collaboration of the kaon S-parameter 
using lattice QCD with Wilson type quarks. An effort is made to control all sys- 
tematic errors except for the quenched approximation. The preliminary result for the 
renormalization group invariant parameter is Bj^ = 0.834(37), which translates to 
B^S(2 GeV) = 0.604(27) in the MS scheme with anticommuting 75. 

1. Introduction 

The kaon i?-parameter remains an important ingredient in current analyses of CP 
violation in the Standard Model It is is defined in QCD with up, down and 
strange quarks by the matrix element of a four-quark operator between kaon states, 



While this matrix element is well-defined, its relation to the full amplitude in the 
Standard Model relies on the hypothesis that the charm quark can be treated as 
a heavy particle. While the quality of this approximation is difficult to assess, it is 
typically assumed to be valid at the five percent level, which sets the scale for the 
precision to be attained by a lattice determination of Bk- 

2. Bk and Wilson type quarks 

Despite recent progress with other fermion formulations, lattice QCD with Wil- 
son type quarks remains attractive because it is computationally cheap and does 
not suffer from mixing between flavour and spin degrees of freedom (unlike stag- 
gered fermions). On the other hand, all axial symmetries are explicitly broken, 
which leads to mixing of operators with opposite chirality, and the possible oc- 
curence of unphysical fermion zero modes in the quenched approximation. In order 
to avoid the latter problem in typical quenched simulations, the masses of pseu- 
doscalar mesons are typically heavier than the physical kaon mass. Concerning the 
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calculation of Bk, the mixing problem is particularly annoying. Decomposing the 
operator 0(v-a)(v-A) = Oyv+AA ^ Ova+av the relevant piece for K'^-K'^ mixing 
is the parity even component, which renormalizes as follows 121 
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(2) 



While the mixing problem can be solved by imposing axial continuum Ward identi- 
ties El it represents a major obstacle for precision results. Note that theparity-odd 
operator component Ova+av does indeed renormalize multiplicatively This can 
be exploited for the computation of Bk as will be explained shortly. 



3. QCD with chirally rotated mass terms 

Axial and vector symmetries can be distinguished according to whether or not the 
quark mass term is left invariant. Hence, a non-standard form of the quark mass 
term also modifies the form of the symmetry transformations. Let us consider the 
continuum theory for a light quark doublet ijj and the s quark including a chirally 
twisted mass term, 

Cf = i^{lp+m + ifiq-/5T^)ilj + s{lp+ms)s. (3) 
After a chiral rotation of the doublet fields 

■0' = exp ^ia75 V', ^' = i/^exp ^1075 , (4) 
and with a chosen such that tana ~ /iq/m, the Lagrangian reads 

C'f ^ii' {Ip+ m') ip' + s{Ip+ rris) s, m' = ^JmFT^. (5) 
The field rotation re-defines the symmetries and maps composite fields, e.g. 

OW+AA = cos(a)Ovv+AA - i sin(Q!)OvA+Av = -»Ova+av (« = 7i"/2). (6) 

In particular, the operators are mapped to each other at "maximal twist" a = 
tt/2 -i^ m — 0, where the quark mass is determined entirely by the chirally twisted 
mass parameter fiq. Hence, by using Wilson quarks with a maximally twisted light 
quark doublet and a standard s quark, the complicated operator mixing problem can 
be by-passed An additional benefit consists in the elimination of unphysical zero 
modes by the twisted mass term. This allows for numerical simulations to get close 
to the physical situation: light, mass degenerate u, d quarks and a heavier s quark. 
However, most results in lattice QCD are currently obtained in the limit where the 
kaon is made out of mass-degenerate d and s quarks. In this way a quenched artefact 
is avoided El but it implies that the zero mode problem is back for the standard 
Wilson s quark. In order to cope with mass-degenerate quarks we also consider a 
different set-up, where the roles of s and u quarks are interchanged. The chirally 
twisted doublet is now ^ — (s, d) and after the chiral rotation one finds 

^W+AA = cos(2q;)Ovv+aa - «sin(2Q;)0vA+AV = -«OvA+AV (a = 7r/4). (7) 
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Again the operators are mapped to each other provided a = 7r/4 <^ m = /iq. For ob- 
vious reasons we will refer to the two set-ups as tt/2 and 7r/4 scenarios, respectively. 
The latter has the advantage that both d and s quark masses can be decreased 
simultaneously without encountering unphysical zero modes. Setting a to tt/2 or 
IT /A requires some parameter tuning using known results for finite renormalization 
constants. In addition one needs the multiplicative operator renormalization con- 
stant. The renormalization problem has been solved non-perturbatively inl^, using 
a finite volume scheme based on the Schrodinger functional and recursive finite size 
techniques 

4. Numerical simulations 

The numerical simulations have been performed using the 0(a) improved Wilson 
quark action and the Wilson's plaquette action for the gauge fields. 4-5 /3-values 
have been chosen in the interval [6.0 — 6.45], corresponding to lattice spacings a = 
0.05—0.1 fm, if the scale is set by rg = 0.5 fm. The lattice volumes range form 16'^ x48 
to 32^x72. Quark masses are tuned to achieve a = 7r/2ora = 7r/4 and pseudoscalar 
masses around or above niK- The analysis of excited states determined safe plateaux 
regions where Bk could be extracted from suitable ratios of correlation functions. 
Finite volume effects were checked at the coarsest lattice spacing {/3 — 6.0), and 
found to be below the statistical errors. In the tt/2 scenario chiral extrapolations 
were performed linearly to the physical kaon mass. In the tt/A scenario the kaon 
mass could be reached by interpolation except at /? = 6.45 where finite volume 
effects at the physical kaon mass would have been non-negligible. The resulting 
values for the Bx are given in figure 1 as a function of a/rQ. Discarding the values 
at the coarsest lattice the data seem to scale very well. A continuum extrapolation 
of 7r/4 data linear in a seems adequate. 

5. Conclusions 

We have performed a benchmark calculation of Bk where, apart from the quenched 
approximation, all systematic errors are under control. These include cutoff effects, 
finite volume effects and the contamination by excited states. A non-perturbative 
renormalization procedure has been employed, and the continuum limit has been 
taken. The preliminary result 

-Bi^ = 0.834(37) ^ sf^(2GeV) = 0.604(27), (8) 

has the desired precision. Within errors it is compatible with results obtained by 
other groups using different lattice regularisations (see ^ for a recent review and 
further references). Hence, further progress will require the inclusion of dynamical 
quark flavours. The recent progress with simulation algorithms for Wilson-type 
quarks ^ will soon allow to compute Bk with a similar precision including the 
effect of light dynamical up and down quarks. 
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Fig. 1. Bk data for both ir/A and tc/2 scenarios. The tt/2 data is sUghtly shifted to the right to 
enhance readibility. Also shown is the Unear continuum extrapolation of 7r/4 data excluding the 
data on the coarsest lattice. 
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